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Q I Abstract: In preparation to the experimental results which will be available in the 

' future, we consider geo-neutrino production in greater detail than in [F. Mantovani et ai, 

arXiv:hep-ph/0309013] , putting the basis for a more refined model. We study geo-neutrino 
■ production for different models of matter circulation and composition in the mantle. By 

^ . using global mass balance for the Bulk Silicate Earth, the predicted flux contribution 

from distant sources in the crust and in the mantle is fixed within ±15% (full range). A 
detailed geological and geochemical investigation of the region near the detector has to be 
performed, for reducing the flux uncertainty from fluctuations of the local abundances to 
the level of the global geochemical error. A five-kton detector operating over four years at 
a site relatively far from nuclear power plants can measure the geo-neutrino signal with 
5% accuracy (Icr). It will provide a crucial test of the Bulk Silicate Earth and a direct 
estimate of the radiogenic contribution to terrestrial heat. 
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1. Introduction 

KamLAND, the Low Energy Anti-Neutrino Detector operating at the Kamioka mine in 
Japan, has demonstrated the possibihty ||T| of detecting geo-neutrinos, the antineutrinos 
from decay chains of radiogenic nuchdes inside Earth. A new window on Earth's interior 
is thus being opened, see, e.g., Refs. |jl], |3|, ||. 

Recently, a reference model of geo- neutrino fluxes has been presented in Ref . . The 
model, hereafter indicated as REF, is based on a detailed description of Earth's crust 
and mantle, and takes into account available information on the abundances of Uranium, 
Thorium and Potassium — the most important heat and neutrino sources — inside Earth's 
layers. 

This model has to be intended as a starting point, providing first estimates of expected 
events at several locations on the globe, see Fig. || ^. In preparation to the experimental 
results which will be available in the future, from KamLAND as well as from other detectors 
which are in preparation, it is useful to consider geo-neutrino production in greater detail, 
putting the basis for a more refined model. In this respect, let us discuss critically the key 
ingredients of REF. 




Figure 1: Predicted geo-neutrino events form U and Th decay chains from i/^ + p ^ e+ + n 
reactions, normalized to 10'^^ proton yr. 

Concerning the crust, the 2° x 2° model of Ref. |^ was adopted for building REF; world- 
average abundances a of radiogenic elements have been estimated separately for oceans, 
the continental crust (subdivided into upper, middle and lower sub-layers), sediments, and 
oceanic crust. Although this treatment looks rather detailed on the globe scale, the typical 
linear dimension of each tile is of order 200 km, so that any information on a smaller scale 
is essentially lost. Since one expects that about one third of the signal is contributed from 
sources within 100 Km from the detector, a better description of the surrounding region is 
needed for a refined estimate. 

^An interactive version of the geo-neutrino event map shown in Fig. ^ wiU be available at the site 
http://www.fe.infn.it/~mantovani/geoneutrini/index.htm . 
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Sandwiched between Earth's crust and metalhc core, the mantle is a 2900 km layer 
of pressurized rock at high temperature. As reviewed by Hofmann in Refs. |^], mantle 
models can be divided into two broad classes, essentially corresponding to the presently 
contradictory geochemical and geophysical evidence of Earth's interior. 

Geochemists have long insisted on a two-layer model, in which the mantle consists 
of a relatively primitive layer below a depth of about 670 km and an upper layer that 
is highly depleted of heat producing elements (panel "a" in Fig. Q). The two layers are 
viewed as separate sources of the Mid-Ocean- Ridge Basalts (MORE), which come from 
mantle regions that have been already depleted in incompatible elements by extraction of 
the continental crust, and of Ocean Island Basalts (OIB), which form by melting of deeper, 
less depleted or even enriched mantle sectors. Also, a more primitive deep layer is needed 
from global constraints, otherwise the amount of radiogenic elements inside Earth is much 
too small with respect to that estimated within the Bulk Silicate Earth (BSE) paradigm 
(see below). 

On the other hand, over the past several years seismic tomography has provided in- 
creasingly detailed images of apparently cold slab descending into the deep mantle, below 
the 670-km boundary. If cold slabs descend into the deep mantle, there must be a cor- 
responding upward flow of deep-mantle material to shallow levels (panel "b" in Fig. |2|). 
If this circulation reaches the bottom of the mantle (whole mantle convection), it would 
destroy any compositional layering below the crust in a few hundred million years (at a 
typical speed of 3 cm yr~^ it takes about 10*^ yr to move down to 2900 km). 

In REF, the Preliminary Refer- 
ence Earth Model Q was used for ^ 
the mantle density profile, dividing 
Earth's interior into several spheri- 
cally symmetrical shells correspond- 
ing to seismic discontinuities. Con- 
cerning its composition, a two-layer 
stratified model was used: for present 
day upper mantle, considered as the 
source of MORB, mass abundances 
of 6.5 and 17.3 ppb for Uranium and 
Thorium, respectively, and 78 ppm 
for Potassium were assumed down to 
a depth Hq = 670 km. These abun- 
dances were obtained by averaging the 
results of Refs 
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Figure 2: Models of mantle circulation, adapted from 
Ref. js): (a) is the traditional two-layer model with de- 
marcation at 670 km and nearly complete isolation be- 
tween upper and lower layers; (b) is a fully mixed model. 



and IJ. Abun- 
dances below were inferred by requiring that the BSE constraint is globally satisfied, ob- 
taining 13.2 and 52 ppb for U and Th, 160 ppm for K. In view of the present debate about 
mantle circulation and composition, it is desirable to have a more general treatment, which 
encompasses both geochemical and geophysical supported models. 
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The Bulk Silicate Earth framework provides a compositional description of the primi- 
tive mantle, the outer part of Earth after core separation and before crust differentiation. 
In this frame, the amount of heat /neutrino generating material inside Earth is determined 
through the following steps: 

(a) From the compositional study of selected samples emerging from the mantle, after 
correcting for the effects of partial melting, one establishes the absolute primitive abun- 
dances in major elements with refractory and lithophile character, i.e., elements with high 
condensation temperature (so that they do not escape in the processes leading to Earth 
formation) and which do not enter the metallic core. In this way primitive absolute abun- 
dances of elements such as Al, Ca and Ti are determined, a factor about 2.75 times CI 
chondritic abundances. 

(b) It is believed, and supported by studies of mantle samples, that refractory lithophile 
elements inside Earth are in the same proportion as in chondritic meteorites. In this way, 
primitive abundances of Th and U can be derived by rescaling the chondritic values ^. 

(c) Potassium, being a moderately volatile elements, could have escaped in the plan- 
etesimal formation phase. Its absolute abundance is best derived from the practically 
constant mass ratio with respect to Uranium observed in crustal and mantle derived rocks. 

The BSE estimates for the total amounts of Uranium, Thorium and Potassium from 



different authors |12, 13, 14, ^ are quite concordant within 10%, the central values being 
w-BSE = 0.8 X 10^^ kg for Uranium, 3.1 x 10^'' kg for Thorium, and 0.9 x 10^^ kg for 
Potassium. These values can be taken — within their uncertainties — as representative of 
the composition of the present crust plus mantle system ^. 

Different models can provide different distributions between crust and mantle, however 
for each element the sum of the masses is fixed by the BSE constraint. This clearly provides 
constraints on the geo-neutrino flux ^ which are grounded on sound geochemical arguments 
or, alternatively, geo-neutrino detection can provide a test of an important geochemical 
paradigm. 

We remark that BSE predictions for the radiogenic contributions to Earth's flow are 
H{V) = 7.6 TW , H{Th) = 8.5 TW and H{K) = 1.8 TW, their sum being about one 
half of the observed heat flow from Earth, He ~ 40 TW. Geo-neutrino detection can thus 
provide a direct insight on the main source of Earth's energetics, see Q. 

The reference model was built so as to satisfy the BSE constraint. However, when 
discussing its uncertainties a wide range of models was considered, from a minimal model 
providing a radiogenic contribution of just 9 TW up to a fully radiogenic model. These 
extreme models, although not excluded from direct observational data, are well outside 
the BSE geochemical constraint. If the BSE global constraint is imposed and only the 



^Note that Uranium absolute abundances in chondritic meteorites are variable within a factor two, 
however the ratio to Aluminum is stable within 10%. 

^Geochemists generally agree on the absence of these lithophile elements in Earth's core, see however 
Refs. 1^, for a different point of view. 

*We are always speaking of the flux integrated over any direction, see Q for a precise definition of the 
geo-neutrino flux. 



-4- 



International Workshop on Astroparticle and High Energy Physics 



Gianni Fiorentini 



crust contribution is left free within the observational uncertainties, then the geo-neutrino 
prediction becomes much more precise. 

Briefly, in this paper we shall address three questions. 

(i) How sensitive are the predicted geo-neutrino fluxes to uncertainties about the mech- 
anism of mantle circulation? 

(ii) Is it possible to test the Bulk Silicate Earth model with neutrinos? 

(iii) Which accuracy is needed concerning the local/regional predictions and the geo- 
neutrino detection? 

We shall restrict the discussion to geo- neutrinos from Uranium progeny, which are more 
easily detectable due to the higher energy. Extension to the other chains is immediate. 

2. Geochemistry, geophysics and geo-neutrinos 

As mentioned in the introduction, the composition and circulation inside Earth's mantle is 
the subject of a strong and so far unresolved debate between geochemists and geophysicists. 
Geochemistry supports the existence of two compositionally distinct reservoirs in the man- 
tle, the borders between them being usually placed at depth near Hq = 670 km, whereas 
geophysics presents evidence of matter circulation extending well beyond this depth. If this 
circulation involves the whole mantle, it would have destroyed any pre-existing layering, 
in conflict with geochemical evidence. In this section we look at the implications of this 
controversy on the predicted geo-neutrino fluxes. 

One can build a wide class of models, including the extreme geochemical and geophys- 
ical models, in terms of just one free parameter, the depth h marking the borders between 
the two hypothetical reservoirs: 

(i) We assume that Uranium abundance in the uppermost part of the mantle is nearly 
chondritic ^ {uu = 6.5 ppb) down to an unspecified depth h. 

(ii) Below h we determine abundances (a;) by requiring mass balance for the whole 
Earth. 

This means that Uranium mass below the critical depth, m>/j, is obtained by subtract- 
ing from the total BSE estimated mass (mesE) the quantity observationally determined in 
the crust (rUc) and that contained in the mantle above h [m^h): 



The abundance in the lower part is then calculated as the ratio of m>/i to Earth's mass 
below h (M^h): 



obtained for h = 25 km {i.e., just below the crust) so that the strongly impoverished mantle 

^We are aware that current estimates of U in depleted upper mantle after crust extraction are in the 
range of 2 to 3.9 ppb Given the uncertainty on these values we prefer to perform calculations using 

the well-constrained CI value. As shown below, the assumption of lower U abundance for the uppermost 
depleted mantle has limited effects on geo-neutrino flux predictions. 
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has a negligible thickness. The traditional geochemical model corresponds to h = h^. As h 
increases, the depleted region extends deeper inside the Earth and — due to mass balance — 
the innermost part of the mantle becomes richer and closer in composition to the primitive 
mantle. 




Figure 3: Generic two- reservoir mantle 
model. Uranium abundance in the upper 
part is fixed at a„ — 6.5 ppb, the critical 
depth /i is a free parameter, abundance in 
the lower part a/ is determined for a fixed 
total Uranium mass in the mantle = 
0.45 X 10" kg. 

Let us discuss in detail a few cases, remembering that the BSE estimate for Uranium 
in the whole Earth is m^sE = 0.8 x 10^^ kg and that the best estimate for the amount in 
the crust is rric = 0.35 x 10^'' kg Q so that Uranium in the mantle is mm = 0.45 x 10^'' kg. 

(a) In the fully mixed model, this quantity has to be distributed over the mantle mass 
Mm = 4.0 X 10^'* kg, which yields a uniform mantle abundance a = 11.25 ppb. We shall 
refer to this model as MIX. 

(b) If we keep the estimated abundance in the uppermost part (a^ = 6.5 ppb) down 
to ho one has the REE model 0. 

(c) Among all possible models, the case h = 1630 km is particularly interesting. Below 
this depth the resulting Uranium abundance is 20 ppb, corresponding to the BSE estimate. 
The innermost part of the mantle is thus primitive in its trace element composition and 
the crust enrichment is obtained at expenses of the mantle content above h. We shall refer 
to this model as PRIM. 

Concerning geo-neutrino fluxes from the mantle, all these models have the same amount 
of heat /neutrino sources and only the geometrical distribution is varied. The largest flux 
corresponds to the model with more sources close to the surface, i.e. to the MIX model. On 
the other hand, the minimal prediction is obtained when the sources are most concentrated 
at larger depth, which corresponds to the PRIM case. 



Figure 4: Uranium abundance in the lower 
part of the mantle ai as function of the crit- 
ical depth h from Earth's surface. 
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From Table ||, the difference between 
the extreme cases is 8%, model REF be- 
ing in between. The abundance in the 
upper reservoir a„ can also be treated as 
a free parameter. If we use an extremely 
low value = 2 ppb |17] down to about 
1200 km and primitive abundance be- 
low, we obtain the minimal prediction 
0.86 X 10^ cm-2 s-^ 

We conclude this section with the 
following remarks: 



Model 


Critical depth h 


Flux 




[km] 


[lO'' cm-2 g-i] 


MIX 


25 


1.00 


REF 


670 


0.95 


PRIM 


1630 


0.92 



Table 1: Mantle contribution to the produced Ura- 
nium geo-neutrino flux. The same Uranium mass in 
the mantle to™ — 0.45 x 10^^ kg and abundance in 
the upper layer = 6.5 ppb are assumed in each 
model. 



uncertainties on the geometrical distribution of trace elements in the mantle can 
change the REF prediction for the mantle by at most ±8%. 

A geo-neutrino detector at a site where the contribution from the mantle is dominant 
(i.e., far from the continental crust) can be sensitive to the mantle compositional 
geometry only if a percent accuracy is reached. 



Since at Kamioka mine or at Gran Sasso the mantle contribution to the total flux 
is about one quarter of the total [Q], uncertainties on the mantle geometry imply an 
estimated error of about 2% on the total flux predicted with REF. 



3. The Bulk Silicate Earth constraint 

So far we have been considering the effect of different geometrical distributions of trace 
elements in the mantle, for fixed amount of these elements within it. Actually the BSE 
model can be exploited to obtain tight constraints on the total flux produced together from 
the crust and the mantle. In fact, with BSE fixing the total amount of trace elements inside 
Earth, geometrical arguments and observational constraints on the crust composition can 
be used in order to find extreme values of the produced fluxes. As an extension of the 
previous section, the maximal (minimal) flux is obtained by placing the sources as close 
(far) as possible to Earth's surface, where the detector is located. 

Concerning Uranium, the range of BSE concentrations reported in the literature is 
between 18 and 23 ppb, corresponding to a total Uranium mass between m(min) = 0.72 
and m(max) = 0.92 in units of 10^^ kg. 

In the same units, we estimate that Uranium mass in the crust is between mc (min) = 
0.30 and mc (min) = 0.41, by taking the lowest (highest) concentration reported in the 
literature for each layer, see Table II of H. The main source of uncertainty is from the 



abundance in the lower crust, estimated at 0.20 ppm in Ref. [19| and at 1.1 ppm in Ref. |20| 



Estimates for the abundance in the upper crust are more concordant, ranging from 2.2 



ppm |21| to 2.8 ppm |12|. We remark that, within this appraoch, the resulting average 



crustal U abundance (occ) is in the range 1.3-1.8 ppm, which encompasses all estimates 
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0.91 ppm, see 



reported in the literature |2|, |T|, |2|, but for that of Ref. |T|], (a, 
Table I ^. 

The highest flux is obtained by as- 
suming the maximal mass in the crust 
and the maximal mass in the mantle, 
mtot(™ax) — mc(max) = 0.51, with a 
uniform distribution inside the mantle, 
corresponding to a = 12.8 ppb. 

On the other hand, the lowest flux 
corresponds to the minimal mass in 
the crust and the minimal mass in the 
mantle, mtot(miii) — mc (min) = 0.42, 
with a distribution in the mantle sim- 
ilar to that of PRIM, i.e. a strongly 
depleted mantle with = 2 ppb down 
to about 1300 km and a primordial composition beneath. 

The predicted fluxes are shown in Table |3| for a few locations of particular interest: 
the Kamioka mine (33°N 85°E) where KamLAND is operational, the Gran Sasso labo- 
ratory (42°N 14°E) where BOREXINO p|] is in preparation, the top of Himalaya (36°N 
137°E), which receives the maximal contribution from the crust, and Hawaii (20°N 156°E), 
a location where the mantle contribution is dominant. 



Reference 


(occ) [ppm] 


Taylor & Mclennan 1985 pl 


0.91 


Weaver & Tarney 1984 |2| 


1.3 


Rudnick & Fountain 1995 [|l9| 


1.42 


Wedepohl 1995 |23] 


1.7 


Shaw et al. 1986 12^ 


1.8 


This work, minimal 


1.3 


This work, reference 


1.54 


This work, maximal 


1.8 



Table 2: Average Uranium abundance in the conti- 
nental crust. 









Himalaya 


Gran Sasso 


Kamioka 


Hawaii 




m 


H 










Crust MIN 


0.30 


2.85 


4.92 


2.84 


2.35 


0.33 


Crust REF 


0.35 


3.35 


5.71 


3.27 


2.73 


0.37 


Crust MAX 


0.41 


3.86 


6.55 


3.74 


3.13 


0.42 


Mantle MIN 


0.42 


3.99 




0.80 






Mantle REF 


0.45 


4.29 




0.95 






Mantle MAX 


0.51 


4.84 




1.14 






Total MIN 


0.72 


6.84 


5.72 


3.64 


3.15 


1.13 


Total REF 


0.80 


7.64 


6.66 


4.22 


3.68 


1.32 


Total MAX 


0.92 


8.70 


7.69 


4.88 


4.27 


1.54 



Table 3: Produced Uranium geo-neutrino fluxes within BSE. Minimal and maximal fluxes 
are shown, together with the Reference values of Ref. Q]. Uranium mass to and heat production 
rate H within each layer are also presented. Units for mass, heat flow and flux are 10^^ kg, TW 
and 10^ cm^^ s^^, respectively. 



At any site the difference between the maximal and the minimal flux predictions are 
of about 30%, the extreme values being within ±15% from the reference model prediction. 



®Note that crust mass and flux ranges found in this paper are tighter than those quoted in Ref. |Q, 
which used the value (acc) ~ 0.91 ppm from Ref. as lower limit. 
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All this shows the power of the BSE constraint. If the total amount of Uranium inside 
Earth is fixed at mesE = (0.8 it 0.1) x 10^'' kg, then the produced geo-neutrino flux at, 
e.g., Kamioka is 

^> = (3.7 ± 0.6) X lO'^cm^^g-i (f^j^ range) (3.1) 

after taking into account the full range of global observational uncertainties on Uranium 
abundances in the crust and uncertainties concerning circulation in the mantle. 

For a comparison, we remind that in Ref. Q, where the Uranium mass was allowed to 
vary in the range (0.4-1.7) x 10^'^ kg, the whole range uncertainty on the produced flux at 
Kamioka was ±1.9 x 10^ cm~^s^^. 

We remark that the error quoted in Eq. ( p.l| ) corresponds to a full range of the predicted 
values. If, following a commonly used rule of thumb, we consider such a range as a zb3o" 
(99.5%) confidence level, we deduce a conventional la estimate: 

^> = (3.7 ± 0.2) X lO'^cm'^g-i (i^j) . (3.2) 



4. The effects of geochemical fluctuations and neutrino oscillations 

The main result of the previous section is that — neglecting regional fluctuations — global 
mass balance provides a precise determination of the produced geo-neutrino fluxes. We 
shall compare this precision with uncertainties resulting from fluctuations of the regional 
geochemical composition, from the available information on neutrino mixing parameters 
and from the detector finite size. 

(a) Actually the Uranium concentration in the region where the detector is located can 
be different from the world average and local fluctuations of this highly mobile element can 
be envisaged. These variations, although negligible for mass balance, can affect significantly 
the flux. In other words, geochemical arguments fix the contribution of distant sources, and 
a more detailed geological and geochemical investigation of the region around the detector 



is needed, the error quoted in Eq. (3.2) providing a benchmark for the accuracy of the local 
evaluation. 

As an example, it has been estimated that about one third (one half) of the geo- 
neutrino signal is generated within a distance of 100 (500) km from Kamioka, essentially in 
the Japanese continental shelf. In REF the world averaged upper crust Uranium concentra- 
tion ttuc = 2.5 ppm was adopted for Japan. In a recent study of the chemical composition 



of Japan upper crust |25| more than hundred samples, corresponding to 37 geologic groups, 
have been analyzed. The composition is weighted with the frequency in the geological map 
and the resulting average abundance is ajap = 2.32 ppm , which implies a 7.2% reduc- 
tion of the flux from Japanese upper crust with respect to the estimate in REF. Larger 
variations occur when rocks are divided according to age or type, see Table and even 
larger differences are found within each group. All this calls for a detailed geochemical and 
geophysical study, with the goal of reducing the effect of regional fluctuations to the level 
of the uncertainty from global geochemical constraints. 

(b) Let us remark that the signal is originated from neutrinos which maintain the 
electron flavour in their trip from source to detector, the effective flux being = ^Pee, 
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\JL \j U.JJ 


Area (%) 


a {ppm) 


Pre-Neogene 


/II '7 


o on 


Pre-Cretaceous 


iU.D 


Z.ll 


Neog-Quat. Igneous rocks 


9/1 1 


9 19 
Z. iZ 


Paleog-Cret. Igneous rocks 


14.1 


3.10 


sedimentary 


39.9 


2.49 


metamorphic 


21.3 


1.72 


igneous 


38.4 


2.48 


Global area weighted average 


99.6 


2.32 



Table 4: Uranium abundance in the upper continental crust of Japan. Groups correspond to 
rock's age or type and quoted abundances for each group are area weighted values, from Ref. 



where $ is the produced flux and Pee = 1 — l/2sin'^(20) is the (distance averaged) survival 
probability. From the analysis of all available solar and reactor neutrino experiments, one 
gets tan^ 9 = 0.41 ± 0.04 at Icr p6|, which implies: 



Pee = 0.59 ± 0.02 (Icr) . (4.1) 

The survival probability is presently know with a 3% accuracy, so that uncertainties on 
the neutrino fate do not mask its sources. 

(c) If Uranium geo-neutrinos are detected by means of inverse beta reaction on free 
hydrogen nuclei {ue + P ^ + n) the event number is 

N = 13.2 6 , , ^ - , (4.2) 

106cm-2s-i 1032 yr ' ^ ^ 

where e is the detection efficiency, Np is the number of free protons in the target and t is 
the measurement time. For a produced flux $ = 4 x 10^ cm~^ s~^ and e = 0.8, one expects 
25 events for an exposure of 10^^ protons yr. Statistical fluctuations will be of order ^/N 
if background can be neglected (this clearly does not hold for Kamioka, due to the many 
nearby nuclear power plants). In order to reach a 5% accuracy — comparable to that of 
the global geochemical estimate — one needs an exposure of 16 x 10^^ protons yr, which 
corresponds to a flve-kton detector operating over four years. 



5. Concluding remarks 

We summarize here the main points of this paper: 

(1) uncertainties on the geometrical distribution of trace elements in the mantle (for 
a fixed mass within it) can change the prediction of the reference model for the geo- 
neutrino flux from mantle by at most ±8% (full range), the extreme values corresponding 
to a fully-mixed and to a two-layer model, with primordial abundance below about 1300 
km. 

(2) By using global mass balance for the Bulk Silicate Earth, the predicted flux contri- 
bution originating from distant sources in the crust and in the mantle is flxed within ±5% 
(Icr) with respect to the reference model. 
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(3) A detailed geological and geochemical investigation of the region within few hun- 
dreds km from the detector has to be performed, for reducing the flux uncertainty from 
fluctuations of the local abundances to the level of the global geochemical error. 

(4) A five-kton detector operating over four years at a site relatively far from nuclear 
power plants can measure the geo-neutrino signal with 5% accuracy. 

(5) This will provide a crucial test of the Bulk Silicate Earth and a direct estimate of 
the radiogenic contribution to terrestrial heat. 
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